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We analyze the scattering and bound state physics of a pair of atoms in a one-dimensional optical 
lattice interacting via a narrow Feshbach resonance. The lattice provides a structured continuum 
allowing for the existence of bound dimer states both below and above the continuum bands, with 
pairs above the continuum stabilized by either repulsive interactions or their center of mass motion. 
Inside the band the Feshbach coupling to a closed channel bound state leads to a Fano resonance 
profile for the transmission, which may be mapped out by RF- or photodissociative spectroscopy. 
We generalize the scattering length concept to the one-dimensional lattice, where a scattering length 
may be defined at both the lower and the upper continuum thresholds. As a function of the applied 
magnetic field the scattering length at either band edge exhibits the usual Feshbach divergence when 
a bound state enters or exits the continuum. Near the scattering length divergences the binding 
energy and wavefunction of the weakly bound dimer state acquires a universal form reminiscent of 
those of free-space Feshbach molecules. We give numerical examples of our analytic results for a 
specific Feshbach resonance, which has been studied experimentally. 

PACS numbers: 03.75.Lm, 34.10,+x, 63.20.Pw, 71.23.An 



I. INTRODUCTION 

An optical lattice is a crystal structure of off-resonant 
laser light, wherein atoms are trapped due to the op- 
tical dipole force. This creates a periodic potential in 
which the motional states of a single atom exhibits many 
of the usual textbook solid state characteristics, such as 
energy band structure, Bloch oscillations and Wannier- 
Stark ladders [TJ [2]. In an interacting many-body sys- 
tem the periodic potential constitutes an ideal setting for 
realizing a wide range of condensed matter phenomena 
such as the superfluid to Mott insulator transition [3j 0] . 
Likewise controlled atomic collisions in an optical lattice 
potential is a promising approach to implement quantum 
gate operations [5]. 

In the context of ultracold atomic collisions a Feshbach 
resonance arises, when a bound state of a closed scatter- 
ing channel is coupled to the scattering continuum in the 
entrance channel [6] . In an applied magnetic field the sep- 
arated atoms are described by eigenstates of the hyper- 
fine and Zeeman-Hamiltonian, and the coupling between 
the different channels is a result of the atomic interaction, 
which is not diagonal in the hyperfme basis. The position 
of the resonance is tunable by varying the strength of the 
applied magnetic field, since the open and closed channels 
have different magnetic moments, and therefore experi- 
ence different Zeeman shifts. When the closed channel 
bound state is tuned to be degenerate with the colliding 
atoms in the entrance channel the scattering cross section 
is resonantly enhanced. A Feshbach resonance may thus 
be used to regulate the atomic interactions, which in the 
ultracold regime is characterized entirely by the s-wave 
scattering length. In addition, since the Feshbach reso- 
nance provides an adiabatic connection between bound 
dimer states and the continuum, a pair of unbound atoms 
may be converted into a molecule by slowly ramping the 



magnetic field across the resonance position. 

To study two-particle physics on a lattice experimen- 
tally, a state with either zero or two atoms at each lattice 
site may be prepared by adiabatic magnetic field sweeps 
across a Feshbach resonance starting from a pure atomic 
cloud [7JIH]. This prepares a state with at most two atoms 
in each lattice well, since sites with more than two atoms 
are eliminated by inelastic collisions. The unbound atoms 
may be subsequently removed by an optical purification 
pulse. The sample thus prepared contains only isolated 
pairs, whose binding can then be controlled by further 
adjustments of the magnetic field. 

A special situation arises in an optical lattice due to the 
band structure. The periodic potential leads to a struc- 
tured continuum with scattering states grouped in energy 
bands separated by band gaps. The depth of the lattice 
controls the width of the bands and their mutual sepa- 
ration. In this paper we concentrate on the lowest band, 
where the band structure then implies that the contin- 
uum has both a lower and an upper edge. This has fun- 
damental implications for both the bound states and the 
scattering states of the system. A counter-intuitive effect 
of the lattice band structure is the existence of bound 
pairs composed of atoms, which repel each other [9} HP]. 
Such repulsively bound pairs exist, when the repulsive in- 
teraction is strong enough to lift the state of two atoms 
on the same lattice site out of the continuum and into the 
band gap, thereby preventing the pair from dissociating 
and making it a bound state of the system [TTJ [I2j [13] . 
In one dimension repulsively bound pairs exist for an 
arbitrarily small repulsion. Near a Feshbach resonance 
further intriguing two-body effects become possible. In 
particular, the Feshbach resonance gives rise to motion- 
ally bound molecules, which are stable only at finite cen- 
ter of mass momentum with respect to the lattice, while 
decaying if at rest. It is thus possible to tune the location 
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of the scattering resonance at a fixed magnetic field by 
controlling the molecular motion [T4] . 

Here we derive analytical results fully characterizing 
both the bound and scattering states of an atom pair in 
an optical lattice near a Feshbach resonance. We repre- 
sent the lattice physics in a discrete model with nearest 
neighbor hopping. To describe the Feshbach resonance 
we adopt a two-channel model with one open channel 
and one closed channel, which is adequate to describe the 
physics of both the scattering continuum and the bound 
states of the system for an isolated resonance [6j [15]. 
Contrary to [14] we include the background interaction 
in the entrance channel, which plays an important role 
both in scattering and in the bound state dynamics. The 
direct scattering process in the open channel interferes 
with the resonant scattering proceeding via the closed 
channel resonance state leading to a Fano profile in the 
transmission probability. In the bound state spectrum 
the presence of an additional bound state leads to avoided 
crossings that are important when considering adiabatic 
transitions induced by ramping the magnetic field. 

The same fundamental physics of a discrete state em- 
bedded in a structured continuum occurs in the Fano- 
Anderson model, which describes autoionizing states in 
atoms [16] and localized electron states in solids [T7] . 
From a many-body perspective the combination of an 
optical lattice with resonant control of the atomic inter- 
actions makes a variety of spin models attainable [18] [19] 
[20]. 

In a previous paper [21] we described the basic results 
of the discrete lattice model of the Feshbach resonance 
physics, and suggested how the salient features of the 
system may be probed in experiments. Here we present 
the full analytic theory of Feshbach scattering in a one- 
dimensional lattice (Sec.[n|). We introduce a generalized, 
one-dimensional scattering length, which is defined at 
both the upper and the lower threshold of the lowest con- 
tinuum band in the lattice. This generalized scattering 
length shows a characteristic Feshbach divergence when 
a bound state of the system crosses either band edge, 
and we identify the magnetic field positions for these di- 
vergences as a function of the lattice depth for a specific 
Feshbach resonance. In Sec. [Til] we consider the bound 
state solutions of the combined lattice-Feshbach prob- 
lem and derive analytic expressions for the binding ener- 
gies and the wavefunctions of the dimer states. We also 
show that the molecules exhibit cooperative tunneling 
with a tunneling rate consistent with that measured in 
experiments [22 , 23 . A spectral analysis of the unbound 
and dimerized atom pair states is presented in Sec. |IV[ 
where we demonstrate that threshold effects lead to a 
non-trivial peak structure of the molecular spectral func- 
tion inside the band. Our discrete lattice model includes 
only the lowest continuum band and neglects couplings 
to higher bands in the lattice. In Sec. [V] we discuss when 
the lowest band approximation is valid and how this re- 
striction may be relaxed. We find that a moderate inter- 
band coupling only modifies our results on a quantitative 



level. Throughout we give quantitative results for a spe- 
cific Feshbach resonance for which our model is valid [24J . 

II. ONE-DIMENSIONAL SCATTERING IN A 
LATTICE 

A. Discrete lattice model 

1. Single particle states in the discrete basis 

We consider atoms with mass m in a cubic optical 
lattice potential V[ at (x) = J2i=i 2 ^o~^R sin 2 (irxi/a) + 

Vq sin 2 (7TX3/a) created by off-resonant laser light with 
wavelength Al- The lattice period is, a = Al/2, 

and the transverse, Vq-, and longitudinal depths, Vq', 
are measured in units of the photon recoil energy, 
Er = h 2 /2m Al- For a cubic lattice the motion separates 
along the three spatial directions, and the Hamiltonian, 
Ho = — {h 2 /2m)V 2 -fViat(x), is diagonalized by products 
of three Bloch functions, (j) nq (xi), one for each spatial di- 
rection i = 1,2,3, where q is the quasi-momentum in 
the range [— 7r/a, 7r/a], and n is the band index. The en- 
ergy eigenstates form the usual band structure, splitting 
the continuum into bands separated by energy gaps that 
decrease with increasing n. The band gaps grow and 
the widths of the bands decrease, as the lattice is made 
deeper by increasing the laser intensity. 

We will assume that the lattice depth along the trans- 
verse directions, Vq~, is much larger than the depth of 

the longitudinal lattice, Vq". This allows us to only in- 
clude the ground state along the two perpendicular di- 
rections, leaving an effectively one-dimensional lattice. 
As the scattering energy is increased consecutive trans- 
verse channels open, however we only consider scattering 
at energies much smaller than the transverse level split- 
ting, thus justifying the neglect of the transverse degrees 
of freedom. A confinement induced resonance occurs at 
the lower edge of the continuum when the (free space) 
scattering length a 3D and the width of the transverse 
confining potential a± = V^^" 1 ) -1 ■> /4 / V are related by 
a ± /a 3B = 1.4603 ... [H [27l M H3 M [31]. However, 
this resonance may be avoided by not choosing the trans- 
verse confinement too tight, and for the quantitative re- 
sults we present here, it plays no role. 

The Bloch waves are delocalized, but can be combined 
into ortho normal localized basis functions, Wannier func- 
tions, which in the longitudinal direction take the form 

^n j (x s )=( — ) / dqe^y nq (x 3 ). (1) 

V27T/ J-vr/a 

In addition to the band index they are labeled by a site 
index, j G Z, specifying at which lattice sites, £3 = ja, 
the Wannier functions are localized. A similar Wannier 
basis can be constructed in the two transverse directions. 
When it is not clear from the context, which Wannier 
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basis we refer to, we will add a superscript _L or || to the 
Wannier functions indicating to which lattice direction 
they belong. 

The dynamics of a single atom in the n'th longi- 
tudinal band, i.e. the kinetic energy and lattice po- 
tential, is characterized by the tunneling amplitude, 
J nm = — (w n j |7Yq I i^nj+m), between lattice sites sepa- 
rated by a distance ma. Here Hq only depends on the 
coordinate along the weak lattice direction, and m is an 
integer. In the transverse directions tunneling can be 
disregarded, due to the much larger lattice depth. 

In this work we consider only the motion in the low- 
est band of the lattice, and we disregard tunneling be- 
tween non-adjacent sites, assuming |Ji m | <C J\\ for all 
\m\ > 1. In the basis of the Wannier states the single 
particle Hamiltonian for motion in the longitudinal di- 
rection then becomes 

h = J(\wij)(u>ij+i\ + ( 2 ) 

3 

with J = Jn. We have set the zero of energy to be the ex- 
pectation value of the energy for an atom in a lowest band 
Wannier state, i.e. (w\j\H\\w\-). The Hilbert space is 
now represented by a discrete basis of the localized Wan- 
nier functions, and accordingly we have translated the 
continuum description of the system into a discrete lat- 
tice model. By our definition E = corresponds to the 
center of the lowest longitudinal Bloch band, and the 
single-particle energies of the discrete lattice model are 
E(q) — —2Jcos(qa). The approximations made are the 
omission of higher bands, and the neglect of tunneling 
beyond nearest neighbor sites. Both of these restrictions 
may be relaxed in a systematic way [32j [33] . 



2. Two-particle states in the discrete basis 

We use the longitudinal Wannier functions as a basis 
for the two-particle states on the lattice. Hence the two- 
particle wavefunction ^(zi, Z2) will denote the amplitude 
of finding the atoms in the Wannier orbitals w m j 1 (#3) and 
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(xs) localized around the discrete positions z\ = j\a 
and Z2 = Considering only the lowest Bloch band 

we have m = n = 1. 

To analyze the motion of the two particles in the lattice 
we introduce center of mass, Z = (^1+^2) /2, and relative, 
z = z\ — Z2, coordinates. For the center of mass motion 
a plane- wave ansatz 



ip(z 1: z 2 )=e lKZ ^ K (z) 



(3) 



is possible, with the first Brillouin zone for the center of 
mass momentum, K, running from —ir/a to ir/a. This 
separation is pertinent, since states with distinct center 
of mass momenta are not coupled by the interactions or 
the lattice, and we can thus obtain separate results for 
the relative motion for each K. Taking the interaction 



potential U to be purely on-site with strength U (de- 
fined in Appendix [A]) , the relative motion of two atoms 
in the open channel is then described by the Hamilto- 
nian, H op = Ho-\-U5 Z fi, where the relative motion lattice 
Hamiltonian, Hq = —2JAf + Ek, contains the discrete 
Laplacian [13] 



Ek 
4 J 



[il>K(z+a)+il> K (z-a)-2tl> K (z)], (4) 



and the center of mass energy, Ek = — 4 J cos (Ka/2). 



B. Entrance channel Green's functions 

The eigenstates of Hq are discrete plane waves, (z\k) = 
y/a/27rexp(ikz). The corresponding energies, ex{k) = 
Ek cos(/ca), are the sum of the single-particle energies of 
atoms with momenta K/2±h, e^(/c) = — 2 J [cos (if a/2 + 
ka) + cos(Ka/2 — ka)). The (retarded) Green's func- 
tion for non-interacting atoms on the lattice is G K (E) = 
[E — Hq + i?7] _1 , where rj is a positive infinitesimal. In 
momentum space it is diagonal 

G° K (E,k,k') = {k\G° K (E)\k>) = 5 E {k _~^ k y (5) 
while its coordinate space form may be found from 



G° K (E,z) 



/ir/a 
-7r / a 



* /a dk 

j a 2ir E — Ek cos(fca) + irj 



(6) 



Due to the delta function interactions we only need to 
consider G K (E,z) = (z\G K (E)\0) . For energies inside 
the band (\E\ < \Ek\) the solution is propagating 



G° K (E,z) 



iex.p(ip\z\) 



(7) 



with pa = cos 1 (E/Ek) 1 while the the solution outside 
the band {\E\ > \Ek\) falls off exponentially: 

G K (E : z) = sgn(E) e ^I^ [-sgn(^)]^. (8) 

Here na = cosh -1 \E/Ek\- For energies above the band 
the sign of G K (E,z) alternates between lattice sites. 

The Green's function related to the open channel 
Hamiltonian, including the interaction, U, is found by 
explicitly solving the Dyson equation 

G K (E) = G° K (E) + G K (E)UG U K (E). (9) 

In coordinate space the result has the simple form 

G K (E,z) 



Gk{e ' z) ~ i-ug k \e,oy 



(10) 



Defining the diagonal elements of the non-interacting 
momentum space Green's function to be Q K (E,k) = 
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[E — ex(k)] 1 the interacting open channel Green's func- 
tion has the momentum space form 



g u K {E,kX) = g° K (E,k,k') 

a Ug° K (E,k)g° K (E,k') 



~2tt 1-UG° k (E,0) 



(11) 



which follows from the matrix elements of Q in the rela- 
tive momentum eigenbasis, using that for an on-site inter- 
action the matrix element, (k\U\k f ) = Ua/2ir, is constant 
for all k'). 



C. Single channel scattering 

With the expressions for the Green's functions the for- 
mal solution of the Schrodinger equation for the relative 
motion becomes explicit. Below we derive the transmis- 
sion probability for open channel collisions and analyze 
the near threshold scattering in terms of a generalized 
scattering length. This section introduces the concepts 
and quantities central to the ensuing discussion of two- 
channel scattering. We first present results relevant to 
scattering of distinguishable particles, and then discuss 
how to treat the case of identical bosons or fermions. 



1. Transmission profile 

Choosing the situation with the incident wave entering 
from the left with quasi-momentum pa = cos~ 1 (E/Ek), 
the scattering state of two particles colliding under the 
influence of the entrance channel potential, J7, is given 
by the solution of the Lippmann-Schwinger equation 

1%(E,z) = e*»+ (z\G%(E)\z')(z'\U\z'')^(E,z''). 

z' ,z" 

(12) 

For our on-site interaction of strength, J7, the scattering 
wavefunction becomes 



$?(E,z) = e ipz + UG%(E, 0)e 



ip\z\ 



(13) 



where we have used that G%(E, z) = G U K {E, 0)e^l z l. We 
note that for the contact potential ip^f(E,z) equals its 
asymptotic form 



f bs (E,K)e 



ip\z\ 



W-oo) (14) 



everywhere, hence we may readily identify the scattering 
amplitude for the background interaction: 



ME,K) 



UG° K (E,0) 
1-UG° K (E,0) 



(15) 



which is unit-less in a one-dimensional system. For the 
one-dimensional problem and distinguishable particles it 



is natural to recast this in terms of transmission and re- 
flection amplitudes. These are defined from the left and 
right asymptotic form of the wavefunction 



tuE^y 



(z -> oo), 



r hg (E,K)e-^ (z - -oo). 



(16) 



By comparing with ( 14 ) it follows that the transmission 
amplitude is tbg = l + /bg> while the reflection amplitude 
coincides with the scattering amplitude, rb g = /b g . The 
resulting transmission probability is 



T bs (E,K) = \t bg (E,K)f = E2 f« u2 E ' E2 - (17) 

As one might expect, the transmission vanishes in the 
limit where \U\ — > oo, and approaches unity as the 
strength of the on-site interaction is diminished. As 
usual, the reflection coefficient is Rb g = l r bg| 2 = 1 — ^b g5 
which is unity at the edges of the band, where propaga- 
tion ceases. The minimum value of the reflection prob- 
ability, U 2 /(E 2 K + U 2 ), is attained in the middle of the 
band, E = 0. We emphasize that the transmission is in- 
dependent of the sign of the interaction. This is a special 
property of the delta-function potential. 

As usual the scattering amplitude may also be related 
to the phase shift, <5b g , of the scattered wave. In one 
dimension the relation is [34] 



f h& {E,K)= 1 -(e 



2i5 hg (E,K) 



1 



)■ 



(18) 



In general there will be two partial waves, corresponding 
to even and odd solutions of the relative motion prob- 
lem. For the delta-function interaction the odd partial 
wave has vanishing phase shift due to symmetry, and the 
reflection and transmission probabilities are then given 
by Rbg = sin 2 (£b g ) and Tb g = cos 2 (<5b g ), respectively. 



2. Identical particles 

The preceding derivation of the one-dimensional scat- 
tering in the lattice assumed that the two colliding atoms 
are distinguishable. If the collision involves identical 
bosons or fermions the scattering state t/j^f must be an 
even or an odd function of the separation z, respec- 
tively, to make the wavefunction either symmetric or 
anti-symmetric under interchange of the two atoms. For 
fermionic atoms in the same internal state this implies 
that the wavefunction vanishes for z = and the contact 
interaction does not produce any scattering. 

For identical bosons the properly symmetrized scatter- 
ing state may be written as 



-ip\z\ 



e 2iS hge ip\z\ ^ 



(19) 



When identical particles collide in one dimension the 
transmission and reflection coefficients have no mean- 
ing, since the incoming and outgoing fluxes are identi- 
cal, however the scattered wave still experiences a phase 
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shift given by (18). Even though the transmission pro- 



file calculated above does not play a role in the scat- 
tering of identical bosons, it may nonetheless be mapped 
out spectroscopically. Consider transitions from a deeper 
bound state i in the molecule to an energy eigenstate 
inside the band induced by applied RF- or photodisso- 
ciation fields, which are described by some transition 
operator T. If the initial state is well localized with 
respect to the lattice spacing, such dissociative transi- 
tions towards the state ( p~9| ) occur with a probability 
\(^(E)\f\i)\ 2 oc \^°i(E,z = 0)| 2 oc cos 2 £ bg . This ex- 
pression, in turn, is given by the transmission probability, 
Tbg(^,i^T), obtained above. 



decaying exponentially with the separation between the 
two atoms and a phase factor exp(±z7r|z|/a), which alter- 
nates between 1 and —1 from one lattice site to the next, 
in accordance with the behavior of the Green's function 
for the relative motion at E > \Ek\ 

As \U\ — » the scattering length diverges, and the 
bound state approaches the edge of the continuum. 
Our intuition from scattering theory in three dimensions 
about the relation between weakly bound states and the 
scattering length thus holds both below and above the 
continuum band. However, since U is proportional to the 
free space background scattering length, <23D, the gener- 
alized one-dimensional scattering length and its three- 
dimensional equivalent are inversely related [26] . 



3. Scattering length 



When the relative quasi-momentum approaches the 
center or the edges of the Brillouin zone, the two-particle 
energies exik) tend to ±\Ek\ and the scattering ampli- 
tude attains a limiting form 



hg(E,K) 



1 



l-ma\E K \/U 



for k — > 0, 



(20) 



with k = k at the bottom of the band and k = ±7r/a — k 
at the top of the band. These limits allow a natural defi- 
nition of a generalized one-dimensional scattering length 
in the lattice 



a bg (K) 



a\E K \ 
U 



(21) 



from the limit /b g —> — [1 + iftftbg] of the scattering 
amplitude as k — > [26 j. An equivalent definition is 
through the derivative of the phase shift with respect to 
the relative momentum 



lim 



d5 hg (k,K) 
8k 



(22) 



in analogy with the usual free space scattering in three 
dimensions. Our lattice scattering length is a generaliza- 
tion of the usual concept, in that it is defined for col- 
lision energies at both the lower and the upper edge of 
the continuum. The sign of a bg at the upper edge of 
the continuum is a convention, which will be justified for 
the two-channel case below. Our generalized scattering 
length depends on the center of mass motion of the pair, 
which is a crucial feature in the lattice, especially when 
we discuss the two-channel case below. 

For U < a bound state is situated below the con- 
tinuum, making a bg positive as expected. In this case 
the pole of (20) lies along the positive imaginary axis, 
k = i\U\/a\EK\, and the analytic continuation of the 
plane wave exp(zfc|z|) is the dying exponential of the 
bound state. Conversely, for a repulsively bound pair 
state lying above the continuum, U > 0, and <2b g is neg- 
ative. The pole of the scattering amplitude is then at 
k = ±7r/a + iU '/cl\Ek\, corresponding to a wave function 



D. Two-channel scattering 

We now consider scattering in the presence of a Fesh- 
bach resonance mixing the entrance channel with an en- 
ergetically closed channel through a coupling, W. This 
gives rise to a set of two coupled equations for the relative 
motion of an atom pair 



K I 



H*\1>$t)+WWg) 



E\il>%). 



(23a) 



(23b) 



Due to the large splitting between the bound states in the 
closed channel the resonance physics is faithfully repre- 
sented by considering only a single bound state, the reso- 
nance state |</>res), of H cl with an energy E res (B, K) with 
respect to the center of the open channel band. This 
amounts to a single-pole approximation for the closed 
channel Green's function 



g£(e,b) 



\4>res}{4>r. 



E-E ies (B,K)' 



(24) 



The resonance state energy E res (B,K) is composed, in 
part, of the free space resonance energy El™ e (B), which 
is linearly tunable with an applied magnetic field, due 
to the difference A/i between the separated atoms and 
the closed channel magnetic moments, and, in part, of 
the displacement of the closed channel relative to the 
open channel Bloch band as illustrated in Fig. [I] At 
all magnetic field strengths we take the zero of energy to 
be the center of the open channel continuum band. In 
free space the bare resonance energy can be related to 
the magnetic field position of the zero energy resonance, 
i?o ree , as outlined in Appendix [bJ 

As a function of continuous position variables 
the longitudinal lattice potential for each atom is 
Vq Er sin 2 (nxsi/ a) with i = 1,2, and therefore the com- 
bined lattice potential for a pair of atoms in the resonance 
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FIG. 1: (Color online) Schematic diagram illustrating the dif- 
ference between two-channel scattering in free space and in an 
optical lattice. The dashed and solid curves are the open 
and closed channel potentials, respectively. For the lower 
pair of potential curves the lattice is absent, and E*l e s e (B) 
is the energy of </> res relative to the open channel thresh- 
old. Applying an optical lattice changes the position of 
the open and closed channel by E a — 2{e\ + 2E^) and 
Em(K) = E™\\K) + 2EF ± (0), respectively, where E\ is the 
average atomic energy and E™ l (K) is the molecular spectrum 
of the first Bloch band both along direction i. Note that the 
position of the open channel thresholds in the lattice depends 
on the quasi-momentum, K, and therefore we choose E — 
to be the center of the open channel band. Hence, the thresh- 
olds are located at ±|£7k|- The difference in energy between 
the closed channel bound state and the center of the open 
channel band in the presence of the lattice is E res (B, K) — 
EllT(B) + E m (K) -E a = E res (B) - 2J m cos(Ka). 



state is 

V^ t (x 3U x 32 ) = 2V^E R sin 2 
+2V<}e r cos 2 



TtX; 



3 ^ 2 / 7t5x 3 
COS I 



a 
a 



sm 



2a 
ttSx 3 
2a 



where we have introduced the center of mass, X 3 = (#31 + 
X32)/2, and relative, Sx 3 = x 3 i — £32, coordinates. Since 
the size of the closed channel bound state is much smaller 
than the lattice spacing, we may consider the limit of 
5x 3 <C a for this state. Consequently, the resonant state 
has a center of mass motion corresponding to a lattice 
potential twice as deep as that governing the motion of 
the individual atoms, while its mass is the sum of the 
two atomic masses. This leads to a tunneling rate for the 
resonant state, J m , which is much reduced in comparison 
with the rate at which the free atoms tunnel through the 
lattice, J m <C J. This bare molecular tunneling leads to 
a modulation of the resonance energy with the center of 
mass momentum of the pair, which we may parametrize 



E res (B, K) = E res (B) - 2J m cos(Ka). 



(26) 



The band averaged resonance energy, E res (B), incorpo- 
rates the relative energy shifts of the atomic and molec- 
ular bands in the lattice. 



1. Transmission profile 

With the explicit expression for the closed channel 
Green's functio n (|24| ) the closed channel part of the scat- 
tering problem (|23b|) may be solved formally: 



\<f>r 



E — E ies (B, K) ' 



(27) 



where we have introduced a coupling matrix element, 
(z, op\W\(j) Yes ) = W5 Zi o, incorporating the structure of 
the resonance state. For a typical resonance in alkali 
collisions the coupling arises due to the difference be- 
tween the triplet and singlet molecular potentials, and 
it is therefore inherently of short range compared with 
the lattice spacing. It is therefore reasonable that the 
coupling only acts between atoms localized at the same 
lattice site. 

The open channel component of the scattering state in 
the presence of the resonance may then be expressed as 



t hg e l 



1 



E-E ies (B,K) 



Z M (E,K) 
(28) 

for all z > 0. Again the relative quasi-momentum 
of the scattering atoms is pa = cos -1 (E/Ek), and 
Z M (E,K) = {faeslWG^^Wlfaes) is the molecular 
self-energy, which for the onsite coupling may be writ- 
ten as: SjiffE, K) = W 2 G^(£,0). 



From the form of the scattered wave (28) the to- 



tal transmission coefficient including both background 
and resonant contributions is found to be of the Fano 
form [16] 



T(E,K)=T hg (E,K) 



(e + q) 2 
e 2 + l 



(29) 



The Fano parameters e = 2(E — E Tes — A)/(HT) and 
q = 2A/ST depend on the resonance shift and width 
functions, which are related to the real and imaginary 
parts of the molecular self-energy, 

T, M (E,K) = A(E,K) -ihT(E,K)/2, (30) 

where the real part 

UW 2 



A(E,K) 



E 2 K + U 2 - E 2 



\E\ < \E K \), (31) 



describes the self consistently determined shift of the res- 
onance position from E res to E* es = £? res + A(i£* es ,iT), 
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and 



KT(E,K) 



E 2 



E 2 

K 



U 2 - E 2 



(\E\ < \E K \). (32) 



is the decay width of the resonance state due to the 
coupling to the continuum. Outside the continuum 
T(E, K) = 0, and the molecules are stable against disso- 
ciation. The real part of the molecular self-energy outside 
the band 



A(E,K) 



sgn(£)>V 2 



E 2 

K 



Usgn(E) 



\E\ > \E K \), 



determines the energies of the bound states for a given 



resonance energy as described in Section III The width 
and shift of the resonance are illustrated in Fig. [2] along 
with the Fano parameters for K = 0. Throughout the 
paper we illustrate our analytical results for a set of pa- 
rameters where U = 1.4 J and W = 2. 2 J. As outlined in 
Appendix |a| this corresponds to 87 Rb atoms in a lattice 
with longitudinal and transverse depths of and 30£?r, 
respectively, in an applied magnetic field tuned close to 
the Feshbach resonance near 414 G 24 . 




c3 
O 



-4 



(b) 










1 




q 



-4 -2 2 4 

E [J] 



FIG. 2: (Color online) The shift, A(E, K) and width, T(E, K) 
of the resonance for K — where the band edges are at 
E = ±4 J (a). The real part of the molecular self-energy is 
continuous across the band edges and has a pole at the un- 
coupled bound state energy, (see section |lll| . The Fano 
parameters, e and q, parametrize the shape of the transmis- 
sion resonance (b). Here E res =2 J. 



of mass motion, due to the weak tunneling of the closed 
channel molecules. This point of total reflection only 
occurs, when the bare resonance state lies inside the con- 
tinuum. Just as in the non-resonant case the reflection 
probability is unity at the band edges. Hence the reflec- 
tion profile has two minima, when \E res \ < \Ek\- When 
the bare resonance state lies outside the band, R(E, K) 
has a single minimum. The two-channel reflection pro- 
files are plotted in Figs. |3|a)-(d) and[4^a) as colormaps 
and in Figs. [3je)-(h) andj^Jc) as solid lines. For the pa- 
rameters used in the quantitative calculations throughout 
this paper J m /J ~ 0.2, and hence the resonance energy 
appears nearly independent of K in Figs. J3^a)-(d). 

The condition for total transmission is eq = 1, corre- 
sponding to E = E res (B, K) — W 2 /U. As E res is ramped 
upwards (downwards) starting below (above) the band, 
this point of destructive interference between the direct 
and resonant scattering amplitudes first appears at the 
magnetic field, where a bound state of the coupled chan- 
nels system enters the continuum to become a scattering 
resonance. As we shall see below, this occurs at two crit- 
ical values where E res = + W 2 /U - exactly the 
points delineating the region, where the condition eq = 1 
can be satisfied. The reflection profile will thus only have 
both a Fano maximum and a Fano minimum if the con- 
ditions |£ res | < \E K \ and |£ res - W 2 /U\ < \E K \ are 
simultaneously met. 

Neglecting the background interaction corresponds to 
the limit Xb g — > 1, or equivalent ly q — > 0. In this case 
the resonance shift vanishes and the reflection probability 
assumes a simple Breit-Wigner-like profile [T4] 



R(E, K) 



h 2 Tl(E,K)/4 



[E - E Tes (B,K)Y 



h 2 rl(E,K) 

4 



(34) 



where hT (E,K) = 2W 2 / v / i^ 
limit of \U\ 
^bg —> and q 



E 2 . The opposite 
oo on the other hand corresponds to 
-> oo, giving rise to a Breit-Wigner-like 
form of the transmission coefficient 



T(E,K) 



h 2 Tl(E,K)/4 



[E — E res (B, K) 



Aoo] 2 



h^{E,K) ' 

(35) 

corresponding to a dip in the reflection probability at 
the resonance, which now occurs at E = E res (B,K) + 
Aoo. The resonance shift, = — W 2 /U, is independent 
of energy and momentum, while the width is given by 
HT^E.K) = 2W 2 V / E 2 K - E 2 /U 2 . Both the shift and 
the width approach zero as the strength of the onsite 
interaction is increased. 



We note that q = —U I \/ E 2 K — E 2 is independent of 
the coupling, W. It is related to the background trans- 
mission probability through the relation Tb g = 1/(1 +<? 2 ), 
which is a consequence of inversion symmetry in the prob- 
lem [25]. 

The reflection maximum of the Fano profile occurs at 
E = E res (B,K) and is nearly independent of the center 



2. Spectroscopic probing of the Fano profile 



The observations regarding indistinguishable particles 



in section II C 2 can be directly carried over to the sit- 
uation with more than one scattering channel. Again, 
the transmission coefficient T(E, K) has no straightfor- 
ward meaning in the scattering of identical bosons. But 
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E[J] 



K [it/ a] 



FIG. 3: (Color online) (a)-(d) Scattering properties and bound states in an extended zone scheme at different values of the 
resonance energy. The colormap inside the band illustrates the reflection probability, R(E,K), ranging from (white) to 1 
(dark blue). The solid lines are the bound state energies of the coupled system, while the dashed curves indicate the diabatic, 
uncoupled bound states with the nearly horizontal line g iving E T es 

and the other corresponding to E® d40|. (e)-(h) Reflection 
coefficient (solid lines) and the background reflection probability, Rh g (dashed lines) for K =0, 0.25, 0.5, 0.75 n/a (from the 
bottom to the top at E = 0). (i)-(l) Closed channel population of the bound states, Eq. (50). Plots are shown for resonance 
energies E res = — 3 J (a,e,i), - J (b,f,j), J (c,g,k), and 5 J (d,h,l). 



following Fano [16] it is clear that the transition proba- 
bility from an initial state i to the stationary scattering 
state ip^? under the action of some transition operator T 
is \T\i)\ 2 oc T(E,K). Therefore, the resonance pro- 
file is accessible spectroscopically by scanning the energy 
of the final state across the continuum. Starting from a 
deeper bound state of the system the dissociated atom 
signal will have the same energy variation as T(E,K). 



3. Scattering length 

The scattering amplitude in the coupled channels case 
consists of a background and a resonant contribution, 
f(E,K) = f hg (E,K) + f res (E,K), where the resonant 
part 



fres(E,K) 



W 2 G%(E,0)/[l-UG%(E,0)} 



E-E ves (B,K)-^ M (E,K) 



(36) 

is obtained by comparing (28) with e ipz (l + /b g + / res ) 
for z — > oo. The phase shiftis related to the scattering 



9 




-5 0_ 5 10 15 

^res [J] 
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FIG. 4: (Color online) (a) Bound state energies (solid lines) 
and reflection coefficient across the continuum band (col- 
ormap) as a function of E res for a fixed K — 0. The hori- 
zontal dashed line is the repulsively bound pair state existing 
above the band for W = 0, while the diagonal dashed curve is 
the band averaged resonance energy, E res . (b) Closed channel 
population of the bound states, (c) Cuts showing the reflec- 
tion probability at E res = 0, -3, 10, and 5J (solid lines from 
top to bottom at E = 0). The dashed curve indicate the 
reflection profile in the absence of the coupling. 



amplitude in the same way as in the single-channel case. 

As in the case of single-channel scattering we define 
the one-dimensional scattering length, af D (K), from the 



*1DJ 



limit of the scattering amplitude, / — > — [1 -\ 
at the extrema of the relative energy, n — > 0. After some 
algebra we arrive at the simple expression 



af D (B,K) = a hg (K) 



1 



W 2 /U 



E res (B,K)-W 2 /UT\E K \\ 

(37) 

for the scattering length, with the upper and lower sign 
corresponding to the top and the bottom edge of the 
continuum, respectively. This lattice scattering length 
is of the standard Feshbach form |35], and reduces to 



the non-resonant expression (21), when the coupling W 



is set to zero. However, the generalized, magnetic field 
dependent scattering length at the top and the bottom 
of the continuum differ, in that they do not diverge at 
the same value of the resonance energy. The motivation 
for the chosen sign convention of the scattering length at 
the upper band edge is guided by the rationale that the 
two scattering lengths then become equal in the limit of 
a deep lattice, where \Ek\ goes to zero. 

It is convenient to define a pair of detunings 
5 ± (B,K) = E res (B,K) - W 2 /U =F \E K \, such that the 
one-dimensional scattering length diverges at the mag- 
netic field strengths, , where S ± (B^K) = 0. We 
thus get two resonance positions, one for scattering near 




FIG. 5: (Color online) (a) Gene raliz ed one-dimensional scat- 
tering length in the lattice, Eq. (37). Solid lines show at 
K = 0, 0.4, and 0.7 tt/cl (from the left to the right), while 
the dashed lines correspond to af B at the same K- values (in- 
creasing from the right to the left), (b) The corresponding 
binding energies below the band (solid li nes) and above the 
band (dashed lines) and the limiting form (46), valid near the 
Feshbach resonances in af B (thin dot-dashed lines). 



the upper band edge, and one for scattering near the 
lower band edge, as shown in Fig. [5] We caution 
that the relative position of Bq and Bq depends on 
the sign of the magnetic moment difference, A/i. If 
A/i > (A/i < 0), then Bq lies above (below) Bq . 
The lattice resonance positions Bq and Bq are shown 

in Fig. ^ as function of the lattice depth, Vq" . They 
are separated by 2|£'x|/A/i and approach each other 

as Vq is increased, since \Ek\ decreases when the lat- 
tice becomes deeper. Furthermore, in the limit where 
both the transverse and the longitudinal lattice depths 
are large, i.e., Vq-,Vq ^> 1, the Wannier orbitals ap- 
proach harmonic oscillator states, and we have that 



W 2 /U A/iA£, while E a 




V^E R +4^fE R tmd 



E m (K) -> Y Vq E r + 2a/Vo £ R . This & ives the follow- 
ing asymptotic behavior of the magnetic field resonance 
position 



< ee +AB[l-/( 2/ )]- 



where £?Q ree is the Feshbach resonance position in the ab- 
sence of the lattice, and the dimensionless function f(y) 
is defined in Appendix [Bj The asymptotic behavior is 
indicated by the black dash-dotted curves in Fig. [6] 

As shown in Sec. IIIB the divergence of the lattice 
scattering length is associated with a bound state enter- 
ing or leaving the continuum. Hence the lattice gives rise 
to a independent shift of the magnetic field threshold for 
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414.4 




momentum. For a given magnetic field, the scattering 
resonance only exists for a range of center of mass motion 
Bloch states, as shown in Fig. |3jc)-(d) [14 j. Outside this 
range there is instead a true bound state of the system. 
Such motionally bound states are stabilized by their ki- 
netic energy, which displaces the pair state into the band 
gap for the relative motion, where the pair cannot disinte- 
grate. Hence, for a fixed magnetic field, where a molecule 
at rest is unstable, accelerating the lattice provides a new 
avenue for tuning the Feshbach resonance. 



FIG. 6: (Color online) Magnetic field position of the scatter- 
ing length divergences as a function of the lattice strength for 
K — 0. Bq and Bq are plotted for = 15 (dashed curves) 
and Vq 1 — 30 (solid curves). In the example A/jl is positive 
so both i?Q~ -curves lie above their respective B^-curves. The 
left axis shows B^ relative to the Feshbach resonance posi- 
tion in the absence of the lattice, and the right axis shows 
the actual magnetic field position for the specific Feshbach 
resonance considered EH. The black dash-dotted curves are 



given by Eq. (38). 



the formation of two-body bound states away from the 
resonance position in free space. This fact, illustrated in 
Fig. [5| was also observed by Orso et al. [36] who solved 
the two-body physics in a one-dimensional optical lattice 
within a single-channel model. The zero-crossing of the 
scattering length occurs for E res = ±\Ek\, which corre- 
sponds to magnetic field values, — W 2 /UA/j j , shifted 
from the resonance positions. 



III. BOUND STATES 

The bound state solutions of the coupled-channels 
problem have energies, E^,, coinciding with the poles of 
the scattering amplitude f(E,K). The pole locations 
satisfy the equation 



W 2 *gn(E b ) 
Eb - E res (B,K) ' 



(39) 



and can be determined as a subset of the roots of a quar- 
tic polynomial. The bound state energies are plotted as 
the solid lines in Figs. |3ja)-(d) andga). If W = the 
entrance channel bound state energy is 



sign(U)JU 2 



E 2 k 



(40) 



with a wavefunction proportional to G° K (Ei ) ,z). This 
single-channel bound state is indicated by the sinusoidal 
dashed line in Fig[3ja)-(d) and the horizontal dashed line 
in Fig.ga). 

An important consequence of the structured contin- 
uum in the lattice potential is that the stability of a 
molecular bound state may depend on its center of mass 



A. Molecular tunneling 

The bound state pairs can tunnel together through the 
lattice with an effective tunneling rate, which depends on 
their binding energy. We now determine the tunneling 
matrix element of the pairs directly from their dispersion 
relation. 

In the absence of atomic tunneling processes, Ek = 0, 
and the energies of the bound states are given by 



Ei 



j=o 



E res (B) + U ± J[E res (B) - U} 2 



W 2 . (41) 



With the inter-channel coupling turned off these zero tun- 
neling bound state energies reduce to the uncoupled pair 
state energies for J = 0; E res and U. With a small non- 
vanishing single-particle tunneling amplitude we may ex- 
pand (39) to lowest order in \Ek/Ei,\ and find 



pJ=0 



El 
2E(= 



(\J\ « \E b \). 



(42) 



Comparing the momentum dependent part of the bound 
state energy with the dispersion relation for non- 
interacting pairs, impair = — 2J pa i r cos(i^a), we see that 
in the limit of a deep lattice the bound atom pairs tun- 
nel through the lattice with an effective tunneling rate 



^Dair — 



2J 2 



(43) 



as one would expect from arguments based on second 
order perturbation theory [371 El EHl EQj E] . We em- 
phasize that the molecule dispersion depends on the ap- 
plied magnetic field. Cooperative tunneling of atom pairs 
consistent with an effective pair hopping matrix element 
given by J pa i r has been demonstrated experimentally by 
measurement of the frequency of population transfer in a 
double well [22] and by demonstrating that the transport 
properties of an attractive Fermi gas in an optical lattice 
depends strongly on the formation of local pairs [23] . 



B. Binding energy 

Instead of the bound state energy it is often instructive 
to consider the binding energy, defined as the distance of 
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the bound state energy from the (upper or lower) band 
edge: 



bind 



= E h 



sgn(E b )\E, 



K\ 



(44) 



The binding energy is positive for a bound state lying 
above the band and negative for a bound pair below the 
continuum. From Fig. [3]^a)-(d) we observe that the bind- 
ing energy increases with increasing \K\ in agreement 
with the single-channel calculation in [36], which only 
discussed bound states with negative binding energy as 
no transverse confinement was included. 

As Eb approaches either of the band edges, the bind- 
ing energy vanishes in a particularly simple way. To see 
this we observe that E 2 — E 2 K w 2 sgn(i^)|i?x|^bind for 
^bind —> 0. Replacing with sgn(Eb)\EK\ in the de- 



nominator on the right hand side in (39) it then follows 
that 
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FIG. 7: (Color online) The effective reduced mass at the up- 
per band edge (solid line) and the lower band edge (dashed 
line) as a function of the center of mass momentum. 



bind 



sgn(E h )U 2 /2\E K \ 



1 



W 2 /U 

E res (B,K)-W 2 /U-sgn(E b )\E K \ 



(45) 



as the bound state energy approaches sgn^^i^l- The 
denominator is exactly the resonance shape of the gen- 
eralized scattering length (37), and the binding energy 



therefore approaches the universal form 



Experimentally the binding energy of the dimer states 
can be measured by modulating the depth of the lat- 
tice. The bound pairs dissociate, when the frequency of 
modulation satisfies a resonance condition, such that the 
dimers can exchange energy with the undulating lattice. 
The resonance is marked by the appearance of unbound 
atoms [9]. 



^bind = 



h 2 



2^(af D ) 2 ' 



(46) 



when the scattering length is tuned to be large with re- 
spect to the lattice spacing a in the vicinity of the Fesh- 
bach resonance. We have introduced the effective reduced 



mass 



h 2 



d 2 6 K 

dk 2 



(47) 



of the atom pair with center of mass quasi-momentum, 
K, in the lattice. Its limits at the top and the bottom 
of the band are =ph 2 /\Ek\cl 2 , with the upper (lower) sign 
applying for k — > ±7r/a (k — > 0). These limits are plot- 
ted in Fig. [7] as a function of K. From the expression 
for the effective reduced mass it is clear that the sign of 
the universal binding energy expression (46) adheres to 



the definition above. The criteria of applicability for the 
quadratic approximation is |^bind| <C minQE'xl, W 2 /U), 
or 1(5*1 < y/2\E K \W 3 /\U\ 5 / 2 . The universal expres- 
sion (|46| is compared with the exact binding energies 
in Fig. ptb) both above and below the continuum. Due 
to the proximity of the zero crossing of a^ D , where (46) 
diverges, the agreement below the band is worse than 
above the band. 

From Eqs. (37) and (46) it is clear that the binding 
energy vanishes when E res = W 2 /U ± \Ek\- This condi- 
tion for the disappearance of the Feshbach bound state 
coincides with the critical resonance energy at which the 
Fano profile develops a point of total transmission, as 
anticipated above. 



C. Bound state wavefunctions 

The formal solution of the bound state problem con- 
sistent with the single-pole approximation for the closed 
channel Green's function is 



Y ( G%(E b ,z)W\<p ies ) 

\4>res) 



(48) 



where the closed channel population, Z5, is determined 



by the normalization condition: (^ p \^ p ) 



WW) 



1. In coordinate space the open channel component of 
the bound state wavefunction therefore becomes 



Z b WG u K {E b ,z) 



(49) 



while the closed channel component of the bound state is 
described by i^f{z) = \fZ\ ) b z $> Here we have used that 
the coupling is of zero range, and that the bare resonance 
state is localized on a single lattice site. The closed chan- 
nel weight of a bound atom pair with energy E^ is 



z b 



w 2 



\E b \ 



[l-UG° K (E b ,0)] 2 (E 2 b 



E 2 K , 



3/2 



(50) 



which is plotted in Figs [3ji)-(l) and|4jb). As the bound 
state approaches the upper or the lower band edge, it 
follows from the limiting form of the binding energy that 
E 2 — E\ w E 2 < (a/af B ) 2 , and consequently Z5 vanishes 
as (U/W) 2 \a/af D \ in the limit where |a/a^ D | — > 0. 
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In momentum space the open channel part of the 
bound state wavefunction is given by 



2tt 1 



G° K (E b ,k) 



UG° K (E b ,z = 0y 



(51) 



which is peaked near the center of the Brillouin zone 
(k = 0) when the bound state energy is negative, and 
peaked near the edges of the Brillouin zone (k = ±ir/a) 
for Ei) > 0, making it possible to discern the bound states 
with Ei) < from those with Eb > in time of flight 
imaging [9]. This reflects that (j)°^{k) is predominantly 
comprised of the fc-states with energies, e^(fc), in prox- 
imity of Eb. For Eb > (Eb < 0) this corresponds to 
the states near the top (bottom) of the band. The height 
of the peaks is reduced as K approaches ±7r/a, where 
all relative motion quasi-momentum states are degener- 
ate in our model. The closed channel component of the 
bound stat e is ind ependent of the relative momentum: 



enters the continuum from above or from below: 



<|*|>^Lmi for|a± D Hoo. (53) 



In this limit the open channel part of the bound state 
wavefunction assumes a universal form 



'ID 



-\ z / a 1D \ 



sgn(E b )] z 



(54) 



The bound state wavefunctions are shown in Fig. [8] in 
both coordinate and momentum space for E res = — J. At 
this resonance energy a bound molecular state exists just 
below the band for K = (lower bound state branch 
in Fig. [3^b)). With a small binding energy, this state 
extends over several lattice sites as shown in Fig [8] For 
the repulsively bound pair state lying above the band 
(upper bound state branch in Fig. |3jb)) the sign of the 
coordinate space wavefunction alternates between lattice 
sites. 
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FIG. 8: (Color online) Open channel components of bound 
state wavefunctions for E re s = — J, where a weakly bound 
molecular state (solid line, narrow bars) at rest is situated 
just below the band. The repulsively bound pair state lying 
above the band (dashed line, thick bars), is peaked at the edge 
of the Brillouin zone, and its sign alternates from one site to 
the next. Th e first column shows the momentum space wave- 
function (51) in units of the second column ip^ p (z). 



The size of the bound state 
Z b W 2 



\A) = 



El 



2[l-UG« K (E h MHE 2 h -Eir 



(52) 



reflects the magnitude of the binding energy of atom 
pairs. In particular, it diverges when the bound state 



D. Sweep experiments 

We have previously shown that a molecule has a fi- 
nite probability of surviving a sweep through the con- 
tinuum band when the applied magnetic field is varied 
across the Feshbach resonance [21 j. The survival proba- 
bility increases the faster the sweep and the weaker the 
coupling to the continuum. Given the existence of mo- 
tionally bound states, an alternative experiment to inves- 
tigate the passage of a molecule through the continuum 
may then be conducted by starting with the applied mag- 
netic field held at a value such that 4 J < \E res — W 2 /U\, 
and molecules at rest prepared in the bound state with 
energy below E Tes (represented by the the lower branch 
of bound states in Figs. [3] and [3]). Accelerating the lattice 
imparts momentum K\ to the molecules. If the magnetic 
field is then ramped adiabatically to a new value, where 
\Eki\ < \E res — W 2 /U\ < 4 J the moving molecules re- 
main stable, while molecules at rest are now subject to 
dissociative decay. With a further acceleration of the lat- 
tice, either in the opposite direction, or in the same direc- 
tion, but letting K cross one or several Brillouin zones, 
the molecules may be swept through the continuum to 
a final state with center of mass momentum Kf, where 
the molecules are stable provided \Ejc f \ < \E res — W 2 /U\. 
The probability of a molecule surviving the center of mass 
momentum ramp through the continuum depends on the 
rate of lattice acceleration and the strength of the cou- 
pling to the continuum states. 



IV. SPECTRAL ANALYSIS 

In this section we will perform a spectral analysis based 
on the full Green's functions of both the open channel 
atoms and closed channel molecules. In general, given 
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a retarded Green's operator, G(E), the corresponding 
spectral function is the imaginary part of the diagonal 
Green's function, A(E,b) = —2Im(b\G(E)\b), in a basis 
| b) of orthonormal states. 



A. The unperturbed case 

For the unperturbed system ([/=W=0) the spectral 
function is 



A° K (E,z) 



2nD° K (E)e(\E K \-\E\), 



(55) 



where the density of states is D° K (E) = |dex(g)/dg| 1 = 
-G° K (E,0)/iir for \E\ < \E K \. The unperturbed system 
has no bound states, and hence the spectral function is 
zero at all energies outside the continuum band. Inside 
the band the scattering states are equally distributed over 
the whole lattice, making the spectral function transla- 
tionally invariant and directly proportional to the density 
of states, which diverges at the band edges, E = 



B. The single-channel case 



from kto—k. These two waves are degenerate and there- 
fore it is natural to expect the average of the two for a 
given energy in the spectral function. We have used that 
the density of states is unchanged for elastic scattering, 
i.e. D Q K (E) is also the density of states for atoms scat- 
tering under the influence of U. 

It is instructive to calculate the total spectral weight 
residing inside the continuum band: 



^A u K (E,z) = l-^-^. (59) 
\E K \ ^ H 



With the addition of the bound state contribution (57) 
the sum rule, J dE A I ^-(E 1 z)/2tt = 1, is then seen to be 
satisfied. 



C. The two-channel case 

When we have two coupled channels there is a spectral 
density for both the open channel atoms and the closed 
channel molecules. We start by addressing the latter. 



Both for the single-channel and the two-channel case 
we will need the entrance channel Green's function, 
G\{E, z, z') = (z\G I ^-(E)\z / ) 1 i.e., the propagator from 
one relative coordinate in the lattice, z' \ to another, z. 
When we consider just a single scattering channel this 
can be obtained from the Dyson equation 



G u K (E;z,z') = G° K (E,z-z') 

, G° K (E,z)UG° K (E, 



l-UG° K (E,0) 



(56) 



where we have used that the unperturbed Green's func- 
tion is translationally invariant. This expression for 
Gk(E, z, z') reduces to Eq. (10) when z or z' is equal 
to zero. 

For the spectral function we have to treat \E\ < \Ek\ 
and \E\ > \Ek\ separately. Outside the continuum 
G^(E, z) is real and we obtain the spectral function 



A u K (E,z) = 2n 



Ue 



-2k°\z\ 



El 



-5(E - E° b ), 



(57) 



where n®a = cosh -1 \E®/Ek\ is the decay constant for 
the bound state. This gives the probability distribution 
for finding a bound state atom pair z lattice sites apart, 
when the bound state energy is E®. 
Inside the band we find 



A%(E,z)=n(\il%(E,z) 



2 )D° K (E). 

(58) 

where ^^(E^z) is the scattering wave defined in Eq. 
(13). Changing the sign of z is equivalent to changing 



the direction of the incoming wave in the scattering wave 



1. The closed channel 

The dressed Green's function for the closed channel is 
found by summing all the terms in the Dyson series based 
on the bare propagator (24). The result is 



G M {E,K) 



1 



E-E ies (B,K)-Z M (E : K)' 



(60) 



where Y. M (E,K) = W 2 G%(E,0) is the molecular self- 
energy pO] ). If E is outside the continuum, Yjm(E,K) 
is real, and the spectral function is only non-zero at the 
poles of Gm, which coincide with those of the scattering 
amplitudes, i.e. the bound states of the coupled channels 
problem ( 39 ) . The spectral weight at the poles is given by 
the corresponding residues, which are just the probability 
of finding the bound state pair in the closed channel, Z^. 



A M (E,K) = 27rZ b 5(E-E b 



(61) 



For energies inside the continuum, we obtain the usual 
structure, 



A M (E,K) 



hT(E, K) 



[E — E Tes (B, K) — A(E, K)Y 



n 2 r 2 (E,K) ' 
4 (62) 

where the energy dependent shift, A, and decay rate, 
T, are defined in Eqs. (31) and (32), respectively. In the 



middle of the band the shift and width are slowly varying 
functions of energy, and the molecular spectral function 
has an approximately Lorentzian shape with a peak po- 
sition given approximately by the resonance energy £? res . 
However, near the band edges significant deviations from 
this simple form arise due to the strong interference with 
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the continuum threshold. This leads to a double peak 
structure of Am, when the bound state is about to en- 
ter the continuum from below, as is evident from Fig. |9j 
which shows Am as a function of the energy and E res for 
a molecule at rest. 



propagator (60), but the residues at the atomic poles 




E re s[J] 



FIG. 9: (Color online) Molecular spectral function as a func- 
tion of the resonance energy for K = 0. Inside the band Am 
is a continuous function given by ( |62[ ). Outside the band the 
delta function peaks at the bound state poles of the propa- 
gator, Gm, are indicated by vertical lines with a height of 
27rZb- Offset below is a contour plot of Am inside the band, 
with bright colors indicating high spectral density. The solid 
lines in this plane are the bound state energies, while 

Eres is 

indicated by the dashed line. 



The molecule lifetime also follows from Fermi's golden 
rule, T(E,K) = 2n\((l) res \W\ilj h ^)\ 2 D (E)/h. Though the 
density of states diverges at the edges of the continuum 
(van Hove singularities) the decay rate goes to zero at 
both thresholds, due to vanishing overlap of the back- 
ground scattering state with the bare resonance state 
when E = ±\E K \. 



2. The open channel 

For the open channel the full Green's function in the 
resonant case reads 



G K {E-z,z f ) = G° K (E,z-z f ) 



G%(E 1 z)^ K (E 1 B)G° K (E 1 -z f \ 



1-E K (E,B)G° K (E,0) 



{63) 



where B) = U + W 2 /[E - E res (B, K)\ is the self- 

energy for the relative motion of an atom pair in the open 
channel. This gives the spectral function 

A°g(E, z) = 2irZ b [WG%(E, z)] 2 S(E - E b ) (64) 

for energies outside the band. The open channel Green's 
functions has exactly the same poles as the molecular 



depend on the separation between the atoms in the pair. 
However, using; the identity E 2 ez e_2/ " N = EG° K (E,0) 
one may show that 



A°g{E, z) = 2tt (1 - Z h ) 5{E - E b ), (65) 

zez 

such that the combined spectral weight of the open chan- 
nel and closed channel at a pole is unity. We must sum 
over all possible separations for the open channel weight, 
since the component of the bound state wavefunction in 
the open channel is distributed over the lattice according 
to (49). 



V. INCLUSION OF HIGHER BLOCH BANDS 

Thus far we have only investigated the dynamics in the 
lowest Bloch band of the lattice potential. In this section 
we analyze the neglected coupling to higher bands. First 
we estimate the size of the second order energy shift from 
the excluded higher bands and derive a criterion for the 
validity of our theoretical model. 



VB 



Next, in Sec. 

discuss how a modest coupling to higher atomic bands is 
expected to modify the central conclusions of this work. 



A. Validity of the lowest band approximation 

In this section we set up a criterion for using the the- 
ory, by comparing the second order energy shift from 
the omitted atomic and molecular exited Bloch bands to 
the relevant energy scales. The dominant contribution 
to the energy shift comes from the bare resonance state 
coupling to atom pairs in higher bands, while the shift 
from open channel atoms coupling to higher molecular 
Bloch bands are several orders of magnitude smaller for 
two reasons. Firstly, the coupling matrix elements be- 
tween an atom pair in the lowest band and an excited 
band bare molecule are much smaller than the coupling 
to a molecule in the lowest band (see appendix |C|) . Sec- 
ondly, the energy splitting between the ground and ex- 
cited molecular bands is larger than the atomic band gap, 
since the molecules experience a deeper lattice potential. 
Hence we approximate the unperturbed wavefunction by 
|0res)> the resonance state in the lowest molecular band. 
For simplicity we will furthermore assume the width of 
the Bloch bands to be small compared to the separation 
from the lowest band, so that the energy can be approx- 
imated by the average energy of the band. 

Under these assumptions the second order energy shift 

is 



E 



|Wl r 



(n,m)^(l,l) 



E n + E„ 



2E, 



(66) 



where the sum is over pairs of atoms in the n'th and ra'th 
Bloch band, Win™ is the inter-band coupling matrix el- 
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ement denned in appendix [Aj and E n is the average en- 
ergy of the n'th band. If the atoms are identical bosons, 
we can take the sum to be over n < m to avoid double 
counting. 

We will require the shift to be small compared to the 
relevant energy scales. If AE^ is comparable to the gap 
between the ground and first excited band 2 (E2 — Ei), 
the higher bands will become populated, and the lowest 
band approximation of our model will break down [42] . 
On the other hand a stronger condition can be derived 
by requiring that the energy shift be small in comparison 
with the width of the ground band, 8J, which sets the 
relevant energy scale for the dynamics in our model. By 
fitting the respective energies from numerical band struc- 
ture calculations at different lattice depths we arrive at 
the condition 



q bg \/A/A/AB\ 
100 a B J \ii B J \1G J 



m 



i) (lOOOnm) < 



QW) 



0.6 ' 



(67) 

where Q(Vq ) = 30(Vq" +4) when comparing 

AE^ w ith 

the band gap, and Q(Vq') = 60 exp(— 0.2Vq ) when re- 
quiring the energy shift to be much smaller than the 
width of the lowest band [43 . The scaling factors are 
the Bohr radius, <2b, the Bohr magneton, //b, Gauss (G), 
the atomic mass unit (amu), and nanometers. If both 
criteria are met, coupling to higher bands can be safely 
neglected. Otherwise, if the energy shift incurred from 
the coupling to atoms in higher Bloch bands is signifi- 
cant with respect to the band width yet negligible with 
regards to the band gap, we might expect quantitative 
changes to our results, but as we discuss below, the qual- 
itative features of our model will remain. In the case 
where the lowest band approximation fails coupling to 
higher bands can be systematically accounted for in the 
tight binding limit [44] [45] [46], leading to an effective 
single-band Hamiltonian with a non-local coupling be- 
tween dressed molecules and unbound atom pairs [T8] . 

Both of the above stated conditions can be fulfilled for 
the Feshbach resonance studied by Syassen et al. in 87 Rb 
near 414 G [24]. In the entrance channel both atoms are 
in the \F = l,m^ = 0) hyperfine state. The resonance 
is narrow, with AB = 18 mG, and has a very small dif- 
ference in magnetic moments between the two channels, 
A/i = 111 kHz/G, which is preferable if Eq. (67) is to 
hold. 

A final condition for our results to be experimen- 
tally relevant is that the variation of the resonance en- 
ergy due to magnetic field noise SB rms is much smaller 
than the bandwidth 8J. With a lattice wavelength of 
830.44 nm [24] the recoil energy is = 3.5 kHz (equiv- 
alent to 0.16 fiK) and for the particular resonance con- 
sidered here E^/Afi =32 mG. For a longitudinal lattice 
depth of Eti we therefore require 5B rms <C 46 mG ac- 
cording to Table [TJ In Ref. [47] the magnetic field noise 
was found to be less than 4 mG. 



B. Corrections due to higher atomic Bloch bands 

As discussed in the previous subsection the major cor- 
rection to our model comes from the coupling to higher 
atomic bands. We now include these couplings into our 
model. For each atomic band we apply the analysis 
of Sec. |IIB| to obtain band specific Green's functions 
G® K n m (E , z) and G^ Cnm (E,z) for the relative motion 
of one atom in Bloch band n and another in band m. As 
for the first band G Q K (E) is the Green's function for non- 
interacting particles and G^(E) is the Green's function 
for particles with an on-site interaction, U. We deter- 
mine G° K ^ n ^ m (E,z) and G^ n>m (£', z) by replacing e K {k) 
by [— 2 J nl cos(Ka/2 + ka) — 2J rnl cos(Ka/2 — ka) + E n + 
E m — 2Ei] in Eq. ([5|, and replacing U by the band 
specific interaction strength U nm in Eq. (10) (for the 
definition of U nm see Appendix [A| . 

We now introduce the coupling W between the first 
molecular band and all the atomic bands. This gives the 
molecular self-energy 



^\WGV{E)W\&, 



2 r U 
lnra| ^K,n,m 



(68) 



which reduces to (|30j) , when we exclude contributions 
from higher bands. Because we only consider energies 
well outside the excited bands, G^ n m is real and slowly 
varying with energy for n, m 7^ 1. Therefore we conclude 
that the lifetime T of the molecules is unchanged and that 
the change in the resonance shift A = Re^M must vary 
slowly. A change of the resonance shift b y an amo unt 
A' will affect some of the results in Sees. 



II D 



III] and 

by 



IV C[ The changes amount to simply replacing E r( 
E res + A 7 , a renormalized resonance energy. 

For E = and K = we obtain for the parameters in 
[24] an approximate relationship for the change in reso- 
nance shift 



A' ~ — 0.004Er,(Vq') -0 ' 1 (1q" l ) 



0.6 



(69) 



due to the higher atomic bands [48 . For the lattice depth 
used throughout this paper A' is only 2% of the band- 
width 8 J, and it is reasonable to neglect higher bands. 

A direct coupling between ground and excited atomic 
bands may be included by renormalizing the local inter- 
actions to include the effect of virtual transitions to ex- 
cited levels through band changing collisions. This leads 
to an energy dependent correction to the coupling pa- 
rameters [33]. But since the correction is inversely pro- 
portional to the gap separating the ground and excited 
bands, it is also expected to be small. 

We have assumed onsite interactions and only nearest- 
neighbor tunneling in all atomic bands which is a ques- 
tionable assumption for the higher bands, but Eq. (69) 



still provides an estimate of the magnitude of the correc- 
tion to the real part of the molecular self-energy. 
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VI. CONCLUDING REMARKS 

We have presented a complete description of the Fesh- 
bach physics of an atom pair in an optical lattice repre- 
sented by a discrete two-channel Hamiltonian. We have 
demonstrated how the structured continuum due to the 
periodic potential modifies the resonance line shape and 
position as well as the properties of the bound states out- 
side the continuum. 

The lattice potential facilitates repulsively bound pair 
states in the gap between the two-particle Bloch bands, 
and near a Feshbach resonance a further class of ex- 
otic bound states may be created, which are only stable 
when their center of mass momentum exceeds a criti- 
cal value, thus enabling us to tune the Feshbach physics 
by controlling the molecular motion. Ignoring the small 
modulation of the resonance energy due to tunneling of 
the bare molecules, the Feshbach molecules are stable 
bound states of the system if their center of mass quasi- 
momentum in the first Brillouin zone satisfies \K\ > K C1 
where 



\E res (B)-W 2 /U\ 
4J 



(70) 



If \K\ < K c the pair state lies inside the continuum 
and decays with a rate T(E,K), which we have de- 
termined analytically. For magnetic fields such that 
\E res (B)—W 2 /U\ > 4 J the Feshbach molecules are stable 
for all K. Note that while the decay rate near threshold 
is proportional to the density of states, which diverges at 
the band edges, T(E, K) goes continuously to zero as E 
approaches the edge of the continuum, since the overlap 
between the bare resonance state and the open channel 
continuum states is suppressed due to the background 
interactions in the open channel. 

An additional effect of the background interactions is 
to shift and modify the resonance profile. With a non- 
zero onsite interaction U the position of the resonance is 
shifted away from the bare resonance energy E res (B, K), 
and the transmission of atoms in the open channel is de- 
scribed by a Fano profile, with a shape determined by the 
energy dependent shift and width of the resonance. The 
energy dependence is most pronounced near the band 
edges, where the coupling to the continuum is strongest. 
This gives rise to threshold effects, when the resonance 
energy is tuned close to the band edge, as evidenced e.g. 
in the double peak structure of the molecule spectral 
function. We have outlined how the energy dependence 
of the scattering can be mapped out spectroscopically by 
inducing dissociative transitions from molecular bound 
states to two-particle scattering states inside the band. 

From the expansion of the scattering amplitude as the 
collision energy approaches the band edges we have de- 
fined a generalized scattering length, which exhibits a 
Feshbach resonance at both the lower and the upper 
threshold of the continuum. The diverging scattering 
length coincides with the appearance/disappearance of 
the Feshbach bound state. In the limit of a deep lattice 



potential the lowest Bloch band narrows, and these two 
resonances coalesce. Since the bare molecules and the en- 
trance channel atoms experience different lattice induced 
energy offsets, the position of both the upper and the 
lower threshold Feshbach resonance is shifted relative to 
the position of the Feshbach resonance in free space. For 
the specific resonance we have considered here, this shift 
exceeds the width of the resonance. 



APPENDIX A: INTERACTION PARAMETERS 

For the interaction between the atoms in the open 
channel we assume that U is a contact potential of 
strength g^ g = 47rft 2 a^g/m, where is the (free space) 
background scattering length. Therefore the on-site in- 
teraction between a pair of atoms in bands n and m be- 
comes 

U nm = .%»• ] { / dXi [w± (Xi)] 4 
i=l,2 ^ 

x j dx 3 [wl (x 3 )wl (x 3 )]\ (Al) 

under the assumption that only the lowest band in the 
transverse lattice is occupied. Since all lattice wells are 
identical, the integrals can be performed over Wannier 
functions centered on any site. 

The coupling to the closed channel is of short range, 
since it is proportional to the difference between the 
atomic triplet and singlet potentials, which have the 
same long range form [6]. In addition, the closed chan- 
nel bound state has a size (typically tens of A), which 
is small compared with the lattice spacing (hundreds of 
nm). We thus assume a purely on-site coupling, which 
in free space would couple molecular and atomic plane 

wave states with a strength g = ^4nh 2 a^AiLiAB/m. 

Here A/a is the difference between the closed and open 
channel magnetic moments and AB is the magnetic field 
width of the Feshbach resonance in the absence of the 
lattice. In the discrete lattice model this results in the 
matrix element 



9 Y[ dx i wf^ ± (x i )w^ (x i ) 2 

i=l,2 J 

x / dx 3 wf ' ll (x 3 )wl (x 3 )w l l (x 3 ) (A2) 



Inm 



for the coupling of a pair of atoms in bands n and m to 
a bare molecule in band I. The matrix elements are zero 
whenever I + n + m is even due to parity of the Wan- 
nier functions. A molecule experiences twice as deep a 
lattice as the atoms, see Eq. (25), and has twice the 



atomic mass. Consequently, separate Wannier functions, 
w nj( x i)i nave t° be calculated for the center of mass mo- 
tion of the resonance state. 

With the exception of Sec. [V] we will consider exclu- 
sively the lowest Bloch band for the atoms and molecules. 
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1 


30 


0.18 


0.043 


1.4 


2.2 


2 


30 


0.14 


0.015 


2.1 


3.0 


3 


30 


0.11 


0.006 


3.1 


4.1 


4 


30 


0.085 


0.003 


4.5 


5.6 


5 


30 


0.066 


0.001 


6.3 


7.6 



TABLE I: Discrete lattice parameters for the 8 Rb Feshbach 
resonance near 414 G. The recoil energy is Er = 3.5 kHz. 



0.478(mC 6 /n 2 ) 1 / 4 , which is a characteristic length scale 
of the long range —Ce/r 6 part of the molecular poten- 
tial [BJ. In the lattice the energy of the closed channel 
resonance state relative to the zero of energy (center of 
the lowest band) is then given by 

E res (B) = AriB-B f ree + ABf(y)}+E m (K)-E a , (B2) 

where the lattice induced energy shifts E a and E m (K) of 
the open and closed channel, respectively, are denned in 
Fig. [J 

APPENDIX C: COUPLING TO HIGHER 
MOLECULAR BANDS 



Therefore, we will only need three parameters to charac- 
terize the dynamics in the lattice: J = Jn, U = Un and 
W = Win- 

As discussed in section |V| our discrete lattice model is 
valid for a narrow Feshbach resonance. As a particular 
example we use the resonance in 87 Rb near 414 G, which 
has been studied by Syassen et al [24]. For this reso- 
nance AB = 18 mG and Afi = 111 kHz/G, while the 
background scattering length is 100.8 Bohr radii. For 
all the plots in this paper we consider a lattice with 
Vq~ = 30 and Vq = 1. This results in interaction pa- 
rameters U/J = 1.4 and W/J = 2.2. Table [j] gives the 
relevant model parameters for this resonance and a range 
of lattice depths. 



APPENDIX B: FREE SPACE RESONANCE 
SHIFT 

In free space the coupling between the collision chan- 
nels induces a shift from the magnetic field B f Y r e e s e , where 
the resonance state is degenerate with the entrance chan- 
nel threshold [E f r r e f (£ r f ^ e ) = 0], to the position B f ree of 
the scattering length divergence. This offset may be ap- 
proximately evaluated by 

~ B l Z e = ABf(y). (Bl) 

Here f(y) — y(l — y)/[l + (1 — y) 2 ] is a dimensionless 
function of the ratio y = / a between the background 
(free space) scattering length and the so-called mean 
scattering length a = (mC 6 /n 2 ) 1 / 4 r(3/4)/2 3 / 2 r(5/4) « 



The molecular Wannier states, Yli=i 2 3 w nlji ( x i)-> 
where rti and 2% refers to the molecular Bloch band index 
and lattice site index along direction z, respectively, con- 
stitute a complete set on which we can expand the atomic 

function /(x) = n^i^I^io^) 2 ] x w L( x sWLo( x s)- If 
we denote the expansion coefficients c^ 1 ^'^ , we have 
Win™ = gcf£ and U nm = g hg \\f\\ 2 , where ||- 1| 2 is the L 2 - 
norm. Now a simple norm argument gives the following 
identity for the neglected molecular couplings to higher 
bands 

s 2 £l«y 2 = ^Hw lnm | 2 , (ci) 

where g 2 /g bg = AfiAB = 0.6E R for the 414 G resonance 
in 87 Rb. For a lattice with V<\ 1 = 1 and = 30 this 
results in the estimate J2i^i |Wm| 2 < 0.0059W 2 u for an 
atom pair in the lowest band (n = m = 1). We therefore 
neglect the coupling of ground band atoms to molecules 
in higher bands in comparison with the matrix elements, 
which couple molecules in the lowest Bloch band with 
excited atom pairs. 
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